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PREFACE ’ . 

The objectives of the contract are to develop a method of using 
LAiTDSAT data tliat vn‘11 ho reliable enough to be used by resource managers 
and planners. This method must be economical, easily-used and readily 
available to the users. Secondly, to set up training sessions to acquaint 
potential u.sers with LAIDSAT data. To find out their needs and involve thsm 

in the continuing research. This involvement will result in the data 

{■ 

collection techniques that will bo accepted bv the potential users. 

The reporting period of this document is from 3/25/76 to G/25/76. 

During this time effort continued to be expended in the development of 
1) computer systems for analysis of the digitally recorded LAflDSAT data for 
forest vegetation , ’"land use and v;atar quality classification and 2) the 
visual analysis of iron mining features and surficial geology. Uork was 
begun for the training of resourcs managers and planners in the application 
of LAflDSAT data. 
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The objectives of this contract are first to develop a inethodology for 
using LAMDSAT data for resource managers or planners and secondly to 
conduct training sessions to acquaint tlie:"? with LAIIDSAT data and appli- 
cations. Three research projects are currently underway , each of v.*hich 
demonstrate an application of LAHOSAT data. These are (1) a classifi- 
cation of forest vegetation and other land uses from digital tapes, (2) 
the compilation of a Quaternary geology map of ilinnesota based on previous 
works and assisted-by visual interpretation of LAIIDSAT bulk imagery and 
(3) projects of regional interest such as visual interpretation of land 
use/1 and cover change for the Mcsabi Range iron mining district and regional 
water quality analysis using digital tapes. 


II 



A. ACCOjlPUSiai.’lTS 

Forest VQnQt.at'ion- and Land Use Cl »issifi cation 

In tirfs fjuorter, all of the planned oroliitn’nary vmrk was conploted - 
including software clevulopnient, training sot selection, and perfonudncG 
of classifications. Although the classification accuracies have not been 
deterifiined yet, general results can be reported at this tiKie. Also, the 
relative value of two algorithms, max inium likelihood and parallelepiped, 
for classifying the vegetation in nortfiern Minnesota can be reported. 

Two new techniques for analyjilng LA!!DSAT data f^-'ere developed and v>ill 
be briefly described. The feasibility of detecting lake Water quality 
iri the study area was considered and the results are included. The value 
of teiiiporally registered data has not been fully evaluated, but some 
initial results i-mre obtained. The abilvty to combine training sets 
for compatability vn'th various classification systems was tested v/ith 
very interesting results, and the implications maybe quite significant. 

The categories of natural resources which 'rill be feasible to map 
from LAilOSAT data have not yet been determined. However, during this 
period classification of the following categories was attempted; 

Water 

Lowland conifer 
Upland conifer 
Mixed forest 
brush and shrub 
brass! and 
Agriculture 
Mined land 
Uncategorized 
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After the accuracy of tiappinij those categories Is deterrn‘no4 sotte clianges 
to the above classification scheiiio may be necessary. 

t fa ter - At this time only one vmter category seems feasible since 
t!m variability In water signatures is large within the lakes, as viell as 
between lakes. This may be due to such factors as turbidity, depth, algal 
growth, tannins, etc., which may have complex interactions. 

2' Lowland Conifers - This category includes black spruce, northern 
white cedar and tamarack. Possibly the tamarack, due to its deciduous 
nature, may be separated if imagery at dates other than May 29 and July 17 
v/arc used. The cedar stands are small and generally occur on the fringe of 
the spruce stands and, hence, from the user's point of view, do not require 
a separate category. 

3. Upland Conifers - The uoland conifer category includes red pine, 
jack pine, and './hite pine. The white pine does not commonly occur in stands 
but as single, typically large, trees scattered throughout the forests in 
Itasca County. Jack pine, red pine and spruce - the latter from the lov/land 
conifer category - are found in relatively pure stands and the signatures 
are obtainable (see Table 1). As a consequence, when signatures for jack 
pine are developed and the classifications are performed, there are cross- 
overs betv/aen jack pine, red pine and spruce. When several jack pine 
training sets were selected it was found that the "overlap" could be pushed 
down - i.e., most of t!ie crossover could be forced betv;een the red pine and 
jack pine, leaving the spruce generally separate. Since the jack pine can- 
not be reliably broken out as a separate category, efforts »;ere made to 
select training sets wwich caused the least overlap with the black spruce 
types. Thus, the upland and lov/land coniferous forest types can bo 
separated but not without some errors. 



The above discussion illustrates a Very important point when consider- 
ing machine processing as a means of interpretinn LA.'IDSAT data. The selec- 
tion of training sots has a tremantloiis effect on the final classification 
results and is subject to bias by the analyst. Also, the importance of 
the investigator's field experience and knowledge of the vegetation types 

r'i 

can significantly affect the classification - the case of jack pins being 
a good example. If tha investigator had not had ground truth for large 
areas in addition to tlis small training area, the first classification 
might have been accepted. The final evaluation, using test sots, would 
undoubtedly show the effects of the spruce-pine crossovers, hut only 
after tlie fact. In this case, adequate ground truth allowed the investi- 
gator to see the problem early in tha process, consider the possible 
alternatives, and adjust the classification system to the liiaitatioris 
set by the spectral signatures. In short, 'why define a category such 
as jack pine if there v/ill be a nigh oercentage of confusion with the 
spectral 1 y-ad j acent categori es? 

■ 4. Mixed Forest - The mixed forest category consists of northern hard- 

wood types (mainly oaks and maples), aspen and birch. These forest types 
were grouped for two reasens. First, tha aspen and birch trees only rarely 
occur as distinct pure stands. In fact, it was nearly impossible to 
locate stands of aspen and birch large enough for use as training sets. 

The northern hardwood stands (although having a distinct signature) are 
not extensive and consequently , v/ere grouped with the aspen/birch (see 
Table 2). 


5. Brush and Shrub - The brush and shrub category includes both up- 
land and lowland types such as alder, hazel , v/illow, leatherleaf and voung 
trees. ^ It includes transition zones between upland and lovM and forest 
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types as v^eU as bot-zoen lakrj shores and the forest* iJhile this catooory 
is genorali includiirj several species » it tnight be adequate for the forest 
manager since it is all nonforested. Hov/ever, the level of detail uould 
hardly be adequate for detailed, wetlaind inventories, although forested 
wetlands (i.e*, spruce bogs) can be mapped. Horo discussion of the 
problems related to mapping wetlands will be given v/hen considering the 
advantages of various algorithms for classification. 

^irassland - .^s with the s!irub/brush category, grassland includes 
both lov/land and upland vegetation. In fact, sedges, cattails and canary 
grass (all wetland species) cosm^i'ise a large* percentage of the area class- 
ified as grassland rather than agriculture. 

7. Agriculture - Agriculture is possibly the most difficult category 
to classify from the LAi!DSAT data. The crops grown in Itasca County in- 
clude mainly alfalfa, timothy, potatoes and sunflowers with the remaining 
land in pasture. The problem results from two factors/ First, the class- 
ification is very dependent upon t!ie date of the imagery. For example, 
freshly plowed fields have considerably different signatures than crops 
later in the season. On a given date, agricultural lands have very dif- 
ferent signatures, but from the forest land manager's perspective, the 
land is all agriculture - bare ground or iriature crop. 

8. Mined Land - Mo evaluation of this category has been dona. 

The vegetation signatures vary with foctors such as crown closure, 

precipitation, pjienology and, perhaps most importantly in tins area, 
intermixing of the plant communities. In Itasca County there are very 
few ’"pure" "stands of any vegetation type, so consi derail a care must go 
into the selection of training sets. If a trainino set is too "pure*', . 
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Table I. Mean Spectral Response from Training Sets of Five Coniferous 
Forest Types, 


Forest Type 


LAUD SAT 

Band 


7^ 

5 

6 

7 

Red pine 

21 .70 

13.91 

3'i.70 

21.30 

Jack pine 

21.67 

Ui.38 

29.63 

17.33 

Black spruce 

22.60 

Wi.OO 

31.'i8 

17.75 

V/hite cedar 

23.50 

16.17 

37.67 

21.33 

Tamarack i 

23.38 

1^.25, 

37 63 

21.13 


Table 2. Mean Spectral Response from Training Sets of Tv.’o Forest Types 
Vlhlch Comprl.se the Mixed Forest Category. 



LANDSAT Band 

Forest Types 

i| 

5 

6 

7 



Aspen/bi rch 

, 

22. Oh 



17.^12 

Northern hardwood 

25.97 



36.15 
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it yill not !« useful in classifyinci that particMlar forest type ar. -‘t 
‘tupicall'' occurs. Furthermore, t!>are arc zones of trarisltian between 
the types which have intermediate sianatures. In areas such as northern 
Minnesota iihere the vetiotatlon patterns are complex, particular Cure must 
also ba givon to selection of the best algorithm for performing the class- 
ifications. Furtharmora, the optimum stratepy used in selecting the train- 
ing set.s is dependent upon which classification algorithm is used. 

Experience in tliis investigation sugfiests that the selection of the 
algorithm depends upon the type of material to i>s classified {i.e., the 
intended application) and the variability in the signatures of the types 
to be classified. Some specific examples will he given later in this 
report to illustrate these points. 

The applicability of the LAi.'DSAT data for detecting patterns within 
v/ater which may relate to water quality was considered. Training sets of 
equal size (each 2% pixels) w.ere selected in open v;ater of five lakes: 

Lake l!innibigoi.hish, Ball Club Lake, Deer Lake, Round Lake, and Jessie 
Lake in Itasca County. These training sets v/are used to ganerate statistics 
(i.e., nean variance, and the range) upon v/hic!i.to classify the antirs 
scene using a level slice algorithm. The means for each of these classes 
are given in Table 3 and are presently being cotiipared and related to ground 
truth, 

l-lhen the level slice algori thm v/as used to perform the classification, 
90.3'd (30,963 of 3'I,2G3) of the pixels classified as water v/ere overlapping 
between the categories defined by the training sets taken from the five 
lakes. This suggests that the variability v/ithin each sanple used as .a 
training set v;as largo, by comnqrison to the variability between the lakes. 
The overlapping. pixels were subsanuently classified usitig a wax imua like- 
lihood algorithm. Thus, in the final classification, all of the water 
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areas were forced into one of tlie five categories of water (s«ni» Table 3). 
The results of tins study shovi a lov.* feasibility to obtain a caafficiOinit 
indicative of the over-all v/ator nuality for an entire lake* If.iagery 
obtained at dates other tfian tliose used in this study raiglit be i.iore useful 
for this application* Various preprocessing transformtions such as band 
ratioing might give different results. P.esearch night also .be done to 
evaluate this application in other parts of Minnesota. 

Although the above conclusions v^are reached, it v/as of intsrast to 
consider the patterns of reflectance v/1th1n lakes. A level slice was 
performed with the Image 100 in each band with various thresholds, iloise 
in the data was very apparent in Band 7 when level slicing tiie Hay 29, 
1973, imagery. The noise was found to be mainly in a very narrow range of 
radiance values, llhen all, of the areas with a radiance value of 3 were 
color coded most of the noise in the data was enhanced. UoisB may ba a 


‘.ignificant factor confounding the classification over the whole scone, 
but it is most apparent in the water. A nearest neighbor or soma type 
of filtering technique might be used before analyzing the data to remove 
the noise, but this was not done in this study. Slicing the same data 
(Band 7, for May 29, 1973) for the range 1 through 2 resulted in alarm- 
ing most of Lake Winnibigoshish, Sand Lake, and iiewstring Lake, but only 
the deepest center part of Ball Club Lake and the shallow edge of Deer 
Lake. When all the pixels \n‘th values of 4 were alarmed none of Lake 
/linnibigoshish v;as colored except near Tauiarack Point. At this level, 
radiance values of 4, the northeast Italf of Ball Club Lake i/as colored 
but none of Bo’wstring Lake. Those oatterns will be comnarGd to nans 
showing lake deptijs. Uhile this has not been comnl cited, the results 


initially suggest that there is no relationship botwGGn'Gither depth 


or turbiditv with the 


Mertaiice values in Band 7 recorded on May 29, 


fi 



Table 3 * Msan Spectral Response from Training Sets of Five Lakes for 
Tv;o Dates. 



• Band 


Band 

9 

Hay 

July 

May 

July 

V/inn i b i cosbi sh ^ 

18.5 

19.6 

9.2 

10.6 

Ball Club 

I8.i| 

19.1 

9^3 ' 

10.3 

i 

Deer 

19.6 

25.6 

9.0 

11. 1 

* 

Round 

■I8.I1 

19.3 

9.8 

10 . 

Jess i e 

17.1 

19.2 

8.9 

10.3 


„g. 






1973. This Is suggested since It Is known that Lake Hinnibigoshlsh Is 
shallow and turbid and Ball Club Lake Is cleaner and deeper. If a 
correlation v/ould have been found between these factors (depth and 
turbidity) and the radiance measured by the tISS, It may be masked by a 
third Interacting factor •• algal growth. 

One application of LANDSAT data for Inventorying »water resources, 
hov/evar, Is unquestionably feasible. The shape, pattern, and extent 
of surface water can be easily mapped. Comparing detailed maps of 
several lakes in Itasca County with the areas classified as v/ater on the 
LAUDSAT data show that some lakes have changed considerably since they 
v/ere mapped. Due to the drought this year, the v/ater level in most lakes 
has dropped. It seems feasible v/ith the LAJJDSAT data to estimate this 
change in v/ater level. 

Two algorithms v/ere used to classify the data Into the nine major 
Cutegories discussed earlier. These include a parallelepiped approach 
and a maximuni likelihood classifier, as implemented on the General 
Electric Image 100 and the Bondix Corporation M-DAS systems, respectively. 
The same data were classified v/ith both systems using the same cateciories 
and nearly identical training sets. It v/as found that the resulting 
classifications v-/ere considerably different when the two systems v/ere 
used. 

When a parallelepiped technique v/as used, less than half of the 
scene v/as classified. The remaining pixels did not fit into any of the 
areas defined In four dimensional feature space, as established by the 
training sets. TIvese results imply that the training sets v/ere ’'atypical" 
examples of the types they v/ere selected to represent. Since all the 
training sets v/ere validated. in the field and carefully located on the 
graymaps, the results v/ere at first nolundGrstand?bTe. Hov/ever, further 
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examination of the photography and more field checking gave a possible 
‘explanation, Pure “classic" stands of each type had been selected for 
use as training sets and were ?iot representative of the types in general. 

The natural’” vegetation patterns of the region are highly variable. Much 
of the total area is covered by mixes of the various types with complex 
transition zones. These patterns result in intermediate signatures which 
v/ere not within the range of radiance values for any of the classes. Using 
the max iwuni likelihood classifier, however, permitted classification of 
more than 9B% of the scene. With this algorithm the pixels of intermediate 
value were forced into one of the categories. 

There are advantages and disadvantages to both classification 
techniques. With the "irst (parallelepiped), it is nearly impossible to 
classify the entire scene and v/ith the second (maximum likelihood), pixels 
are apt to be classified incorrectly. However, there seems to be a logical 
approach to the diTetnma considering the tradeoffs. First, the scene should 
be classified using a parallelepiped technique. The results of this 
cl assif icatiop-shoul d be quite accurate. The signatures classified with 
this technique would hove to be quite similar to those for the "classic" 
training sets and should not include the mixtures and transition zones. 

Once these classifications are performed they could be printed with various 
symbols all in one color. Then the remaining pixels can be classified 
using the maximum likelihood algorithm. The results of this second classi- 
fication cf^uld be printed with the same symbols as used in the parallelepiped 
technique but this time in a different color. The results of such a two 
step procedure v/ould permit mapping nearly the entire scene at two levels 
of confidence (the accuracy could be determined for areas classified with 
each technique individually). The user of such a map would thus know for 
which areas photography and field checking are most needed. 
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Some general statements can be made concerning the relative merits of 
the maximum likelihood and the parallelepiped techniques. Classification 
of large distinct natural resource types can be done using the Image 100 
system. Ildv-zever, when small units of vegetation are to be mapped, or for 
transition zones, the maximum likelihood classifier is far superior. In 
Itasca County, water is very distinct and occurs in large fairly homogeneous 
blocks permitting very adequate mapping with the parallelepiped technique. 
When the same algorithm was used for typing v/etland vegetation, v/hich often 
occurs in narrow bands, the results v/era very poor. In this case pixels 
along the shores of lakes were nearly, fill classified as wetland vegetation. 
Using the maximum likelihood technique the shoreline pixels were usually 
mapped correctly (either as water or the same as the adjacent vegetation 
type - on the land). That the selection of training sets seems to be more 
critical with the parallelepiped technique is understandable since training 
sets including the whole range for each type must be included (i.e.. 
Variable density, age, and understory) when developing. the themes. Thus, 
several training sets may be required in defining the signature for each 
vegetation class while, with the maximum likelihood classifier, only one 
training set representing the center of the range of signatures for each 
type is needed'. V The selection of training sets is the most important step 
in classifying LAM DSAT data by any supervised machine processing technique. 
Human bias in selecting the training sets can significantly influence the 
classification results. !; 

The value of using temporally registered data is difficult to deter- 
mine. Tv/o LANDSAT scenes for Itasca County (May 29, 1973 and July 17, 1974) 
were registered to each other by techniques used at the Jet Propulsion 
Laboratory, Pasadena, California. Training sets were defined and three 

'■ I. 

classifications were performed using the Image 100 system. A portion of 
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the scene was mapped first using only bands 5 and 7 from the Mav imagery, 
then using the same bands on the July imagery. Finally, using bands 5 and 
7 from both dates simultaneously, a third classification was performed. 

The accuracies have not yet been determined but the three color-coded maps 
are quite different although exactly the same training sets were used in 

each case. Thus it can be concluded that the degree of crossovers >”>d types 

\ 

of crossovers change with the phenology of the vegetation. After t‘- 'accur- 
acy for each of the three classifications is determined more comments con- 
cerning this problem may be possible. 

Bendix Corportion recently completed a project with the Bureau of Out- 
door Recreation, U.S. Department of the Interior, for mapping several land 
use types along the Mississippi River and includes all of Itasca County. 
Considerable time was spent evaluating the maps produced. Many of the 
categories seemed ill-defined and very inadequate for forestry applications. 
However, the personnel at the Bureau of Outdoor Recreation found the maps 
to be well suited to their needs. In this project, Bendix developed 35 
original categories of land cover and subsequently grouped the categories 
to form the classes desired by the Bureau of Outdoor Recreation. Since the 
original 35 cata.qories were stored in the computer it was possible to simply 
regroup the categories by changing the color coding. This was tried and 
hard copy prints of these classifications are now being produced. When these 
products are obtained the accuracy will be evaluated and compared to the 
results from the Image 100. 

Tv/o new techniques were developed for analysis of LANDSAT 70mni imagery 
using a density level slicer. _First, a means of determining the distribution 
of the gray levels in a scene ;(i.e. , histogram of radiance) was tested and 
found to be accurate and repeaiable. Second, and unsupervised "quasi- 
clustering" technique for classification using graph paper and a pocket 
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calculator was found to be useful In analysis of the LANDSAT data, A 
manuscript describing both of these techniques In detail Is being prepared, 
and should be finished within the next quarter. 

A. ACCOMPLISHMENTS - cont. 

Regional Development Commission Projects 

During the reporting period the interpretation of all 28 Hesabi Iron 
Range Quadrangles was completed. Final drafting and editing was also com- 
pleted, Twelve 1 : 24 ,000 scale transparent quad overlays are currently being 
evaluated by ARDC personnel. The remaining 16 quadrangles are now In the 
process of photo reproduction. Statistics on the types of land cover changes 
are nearly complete. Total changes from 1969 through 1976 are expected to 
amount to over 18,000 acres, about half of which are expansion of mine 
features. The total changes represent roughly two percent of the 28 quad- 
rangle are In 6 years. 

Most of the problems associated with the use of computer LANDSAT tapes 
for water quality analysis have been over come- and digital tapes covering the 
Douglas county area In Region 4 Regional Development Commission have been 
ordered. 

Quaternary Geology Project 

Progress on the compilation of the Quaternary Geologic map of the State 
of Minnesota has permitted the location and classification of bog sites 
throughout the state. Sufficient data has been collected to proceed v.*1th the 
Identification and location of glacial features. 

The bog areas of the state were generally located with the assistance 
of the Minnesota Soils Atlas. A black and white (positive) print of May 
Imagery at a scale of 1:500,000 from band seven was used to locate the bogs. 
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They were identified visually by darker tones associated with v;et areas. The 
<jrassy bogs reflected more light and did not have the darker tones associated 
with wetness. In these cases the shape and location as well as the apparent 
flatness of the area v/as used to define the bog. The bogs were classified by 
their association to the surficial drainage pattern. 

Technology Transfer Sessions 

Work continues at two levels for the training of resource managers 
and planners in the applications of LANDSAT data. 

Assistance and information was given on a one to one basis to the 
staff of Region 5 Regional Development Commission and Region 7 a Regional 
Development Commission. As a result of these meetings, the Region 5 
Regional Development Commission is writing a contract with the Bendix 
Corporation for the purchase of regional land use maps generated from 
LANDSAT data. •- 

Secondly, work was begun on the text and 35mm slides for a multi- 
media (slide-tape) show introducing the various resource inventory pro- 
grams, including LANDSAT, in Minnesota. 

r- 

B , Problems 

There were no significant problems during this reporting period. 

C. Significant Results 

As mentioned previously in this report, the application of LANDSAT 
data in digital tape format is very feasible. The shape, pattern and extent 
of surface water (eg lakes) can be readily mapped. Comparing detailed maps 
of several lakes in Itasca County with the areas classified as water by the 
LANDSAT data show that some lakes have changed considerably since they were 
lapped. 
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Due to the several droughts this year (1976), the water level in most 
lakes has dropped. At this time it seems feasible that LANDSAT digiWl tape 
data to estimate lake water level change, due to the 1976 drought conditions. 

D. Publications 

There were no publications during this period. 

E. Recorrmendations 

There are no recommendations at this time. 

% **• # l> 41 « «• « - I ^ 

H. Aircraft Data 
None at this time. 

I, Work Planned for next quarter 

Forest Vegetation and land Use Classification Project 
The products of the Image 100 work v/ill be enlarged and the classi- 
fication accuracies will be estimated by sampling several areas and comparing 
them with the ground truth. The Bendix products will also be evaluated. 

Once these statistics are known, the products will be duplicated and de- 
livered to the user-cooperators for their evaluation. Line printer class- 
ification maps will be produced and converted to transparent overlays for 
the areas included in several Ih USGS topographic quadrangles. These material 
will also be given to the field cooperators for their evaluation. Plans have 
been made for working with UfllVAC Corporation to determine if various pre- 
processing operations will enhance the imagery. Detailed descriptions of the 
two new techniques developed for analysis of LAtlDSAT imagery will be written. 
Finally, the feasibility of using a film recorder to classify vegetation types 
in the study area will be evaluated. 





ReglonaV Development Commission Projects 

In the next reporting period the lake classification procedures v-/in 
he developed and attempts vnll be made to apply these procedures to lakes 
in the Douglas County area. Arrowhead Regional Development Coiiraission 
personnel will also complete their evaluation of the Iron Range land cover 
change maps. In addition Department of Natural Resources personnel in the 
nibbing office have expressed strong Interest ir che quads. A set will be 
provided to them for evaluation in their work, by the end of July. 

Quaternary Geology Project 

The next step in the compilation of the Quaternary Geology map will be 
to visually recognize glacial features {moraines, lake plains, outwash plains, 
etc.) on the LANDSAT imagery that were identified in previously published 
books and maps. As the features are recognized on the LANDSAT imagery they 
will be located and delineated on the draft map. 

During the next reporting period it is anticipated that the character- 
istics of LANDSAT imagery which facilitate the recognition of glacial 
features will be understood. This will be the final step in compiling the 
Quaternary Geologic Map of Minnesota. The maps will thgn be prepared for the 
printing process. Later the techniques used in reco'gniziing glacial features 
will be presented to people and agencies who wish to use remote sensing in 
evaluating their resources. • 
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REGIONAL DEVELOPMENT 
COMISSION name 

1 Northwest 

2 Headwaters 

3 Arrowhead 

4 Ifes t Cen tral 

5 Region Five 
6E Six East 

6W Six West 
7E East Central 
7W Central Minnesota 
8 South«e55t 
S Region Nine 

10 Southeastern 

11 Metropolitan Council 
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